Discretisation effects in the topological susceptibility in lattice QCD 
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We study the topological susceptibility, Xi m QCD with two quark flavours using lattice field 
configurations that have been produced with an 0(a)-fmproved clover quark action. We find clear 
evidence for the expected suppression at small quark mass, m q , and examine the variation of \ 
with this mass and the lattice spacing, a. A joint continuum and chiral extrapolation yields good 
agreement with theoretical expectations as a, m q — » 0. A moderate increase in autocorrelation is 
observed on the more chiral ensembles, but within large statistical errors. Finite volume effects are 
negligible for Leutwyler-Smilga parameter sls > 10, and no evidence for a nearby phase transition 
is observed. 



o 
o 

(N 
S3 

a 
>—» 



> 
in 



o 
o 

-f— > . 

i 

<D . 

^3 ■ 



X 



PACS numbers: 11.15.Ha, 12.38.Gc 



I. INTRODUCTION 

QCD is distinguished from the constituent quark mod- 
els of hadronic physics by the presence of light sea quarks 
in the vacuum and one of the principal goals of lattice 
field theory is to provide an accurate theoretical descrip- 
tion of their effects. 

The topological charge, Q, and its associated suscep- 
tibility, Xi are especially sensitive to the presence and 
properties of the sea quarks. Chiral perturbation the- 
ory describes their variation with the quark mass, m q . 
By comparing lattice Monte Carlo measurements with 
this expectation, we may understand how near we are to 
achieving the above aim. 

Initial work on unimproved sea quark formulations 
(Wilson and staggered 0, H, 0, H E ) showed only 
weak evidence for the expected chiral (m q — > 0) sup- 
pression of x relative to the quenched (m q — ► oo) the- 
ory at an equivalent lattice spacing, a (see H, 0, 
for reviews). A clearer chiral reduction was seen for im- 
proved fcrmion formulations, probably as a result of bet- 
ter chiral symmetry at finite lattice spacing [jj. 0(a)- 
improvcd Shcikholcslami-Wohlcrt (SW, or "clover") Wil- 
son fermions showed a strong sea quark effect in a st udy 
at quasi-fixed lattice spacing, a ~ 0.1 fm 0, 0, [n|, 0, 
IT3I Ibi; mean field improved SW quarks also showed an 
effect 15]. Improved staggered fermion studies have also 
shown clear suppression of x over a range of lattice spac- 
ings 0.09 < a < 0.12 fm 0, [11113 Studies using lattice 
sea quark formulations with chiral properties more com- 
parable to those of the continuum have so far been limited 
to two dimensional models 

In the light of the improved staggered results, it is in- 
teresting to examine the effects of finite lattice spacing 
on the topological susceptibility from C(a)-improved SW 
quarks. This is useful not only from the theoretical view- 
point, but also when comparing the relative costs of pro- 
ducing independent field configurations in Monte Carlo 
simulations (the topological modes are expected to be 
amongst the slowest to decorrelate) . 

In this paper we extend the analysis of two flavour 



0(a)-improved SW results at a ~ 0.1 fm (included in 
this work with due acknowledgment) to cover a com- 
parable range of lattice spacings to the improved stag- 
gered fermion study, using lattice QCD ensembles pro- 
duced by the UKQCD and QCDSF collaborations. Mov- 
ing to small quark masses (the lightest yet produced for 
Wilson-like fermions), we find strong evidence for the 
chiral suppression of x- Finite volume effects are shown 
to be under control, and no strong autocorrelation is seen 
for the topological charge. We perform a joint continuum 
chiral and continuum extrapolation, and find good agree- 
ment with theoretical expectations. 

The structure of this paper is as follows: in Section ITT1 
we review the theoretical expectations for the chiral vari- 
ation of the topological susceptibility in the continuum 
and on the lattice. We describe our measurements in 
Section IlIII and discuss their variation with m„, a and 



the lattice volume. 
Section CVI 



We conclude with a discussion in 



II. THE TOPOLOGICAL SUSCEPTIBILITY 

Four-dimensional, Euclidean gauge fields can be clas- 
sified by an integer-valued topological charge, 



Q = J d 4 x^e^ crT F tlu (x)F crT (x) E Z 



(1) 



In the absence of a 0-term, (Q) = and the topological 
susceptibility, x = (Q 2 ) /V, has a finite, non-zero limit 
as the volume V — ► oo. 

As the mass of the sea quarks is reduced the topological 
susceptibility is suppressed, the leading order behaviour 
being [HEIliS 



X(ml) 



(2) 



for two degenerate flavours, with corrections at 0(m 4 ). 
This relation should hold when V is large enough that 
chiral symmetry is not restored, i.e. 



£ls = m q Ty ~ (Um^) 2 V > 1 



(3) 
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TABLE I: The ensembles including the volume L 3 T in lattice units, the number of configurations analysed, their separation in 
HMC trajectories, physical parameters and various measures of the physical volume. 



label 


P 




L 3 T 


A^conf, sepn 


ro 




(r m 


,) 2 


L/r 


Lrrir: 




?i 


5.20 


0.13420 


16 3 32 


514, 10 


4.077 


(31) 


5.67 


(7) 


3.9 


9.3 


150 


U2 


5.20 


0.13500 


16 3 32 


789, 10 


4.754 


(40) 


3.87 


(5) 


3.4 


6.6 


55 


U3 


5.20 


0.13550 


16 3 32 


830, 10 


5.041 


(40) 


2.15 


(3) 


3.2 


4.6 


24 


U4 


5.20 


0.13565 


16 3 32 


280, 10 


5.246 


(51) 


1.67 


(4) 


3.0 


3.9 


16 


U 5 


5.20 


0.13580 


16 3 32 


279, 10 


5.320 


(50) 


1.22 


(4) 


3.0 


3.3 


11 


qe 


5.25 


0.13460 


16 3 32 


194, 10 


4.737 


(21) 


5.44 


(4) 


3.4 


7.9 


79 


11- 


5.25 


0.13520 


16 3 32 


822, 10 


5.138 


(45) 


3.90 


(5) 


3.1 


6.1 


41 


Qs 


5.25 


0.13575 


24 3 48 


91, 10 


5.430 


(60) 


1.99 


(4) 


4.4 


6.2 


85 


Ug 


5.26 


0.13450 


16 3 32 


415, 10 


4.708 


(52) 


5.73 


(9) 


3.4 


8.1 


85 


uio 


5.29 


0.13400 


16 3 32 


397, 10 


4.813 


(45) 


7.70 


Cn) 


3.3 


9.2 


105 


qn 


5.29 


0.13500 


16 3 32 


587, 5 


5.227 


(37) 


4.84 


(5) 


3.1 


6.7 


47 


qu 


5.29 


0.13550 


12 3 32 


834, 5 


5.756 


(33) 


3.53 


(3) 


2.1 


3.9 


10 


qi3 


5.29 


0.13550 


16 3 32 


911, 5 


5.560 


(30) 


3.30 


(3) 


2.9 


5.2 


25 


qu 


5.29 


0.13550 


24 3 48 


405, 5 


5.566 


(20) 


3.30 


(2) 


4.3 


7.8 


127 


qw 


5.40 


0.13500 


24 3 48 


653, 2 


6.088 


(32) 


6.03 


(5) 


3.9 


9.7 


162 



where £ is the chiral condensate [21] 38]. 

As m q — > oo, x approaches the (constant) quenched 
value, X^ qu \ from below. Higher order corrections to 
Eqn. J3| thus introduce a negative curvature at some 
intermediate quark mass. Various interpolating ansatze 
between the chiral and quenched regimes have been pro- 
posed H, 0, |2l|. In this paper we shall focus on one 
particular form, the large- N c form |2l|, as motivated for 
SU(3) by Diirr 



4 

r oX 



c 3 co(r TO 7r ) 2 
C3 + c (r m x ) 5 



with c = (r / T ) 2 /4 , (4) 

C3 



r 4 x (qu) 



Assuming the physical values of ro = 0.5 fm [22| and /„. = 
93 MeV, then Co = 0.0139. This formula has been seen 
to work well in QCD at fixed lattice spacing. 



A. Lattice QCD 

Eqn. (J2J also holds for lattice QCD at finite lattice 
spacing for fermions with sufficiently good chiral proper- 
ties |23|,|24j. In general, however, the lattice topological 
susceptibility is related to the continuum by additive and 
multiplicative renormalisation factors, Z and M: 



r V lat > 

'OA 



^ 'OA 



through the coupling, g 2 ). The upshot of this is that 
discretisation effects act to reduce x below its contin- 
uum limit at large m q , but to increase it in the chiral 
limit. A recent discussion can be found in ,17]. For- 
mally, before examining the sea quark mass dependence, 
we should first perform a continuum extrapolation of x 
at fixed m q to remove the effect of these renormalisation 
factors. This is, however, not possible with the ensembles 
we have available (especially if we do not admit interpo- 
lation in m q of measurements at finite a). To attempt a 
full analysis of our data we must assume that the con- 
tinuum and chiral limits commute sufficiently well that 
we can either perform the m q — > limit first at fixed a, 
or that we can carry out a joint chiral and continuum 
extrapolation using a single formula. 

The first approach was used in H,0|, fitting Eqns. © 
and Q to data at a ~ 0.1 fm. Constraints of simulation 
mean that this study cannot contribute many ensembles 
at comparably light sea quark mass but different a. We 
cannot, therefore, extend the above method by carrying 
out chiral fits at other lattice spacings so we can look 
at discretisation effects in the fit parameters. Instead 
we follow the second approach. Allowing C0.3 to have a 
leading order lattice spacing correction (quadratic for the 
lattice action considered here) 



M . 



(5) 



co(^o) = coo + 



cqi 



C3(m) = C30 + 



C31 



(6) 



These competing factors are in general functions of the 
lattice spacing and the quark mass (partly implicitly 



in Eqn. ijlf. we obtain: 



3 



4 

r X 



C3QCoo(r TO 7r ) 2 

C30 + c on (r a m 7T ) 2 



ki 



(C30 



(roTO„-) 2 



coo(roniir) 2 ) 2 
I 



(coiC3 + c 31 cl (r m n ) 



(7) 



We have added a general constant, &i > 0, to allow x 
to have a non-zero chiral limit at finite lattice spacing, 
as might occur in the absence of exact zero modes in 
the lattice Dirac operator A check of our results is 
consistency of the fitted parameters with the results for 
ro — 5 from 



III. MEASUREMENTS 

We have calculated the topological charge and sus- 
ceptibility using ensembles generated by the UKQCD 
pHil and QCDSF [23 EH collaborations. De- 
tails of these ensembles are given in Table [Q (taken m 
the main from 28] ). The lattice spacings lie in the range 
0-08 < a < 0.12 fm. The SU(3) gauge fields are gov- 
erned by the Wilson plaquette action, with two flavours 
of SW fermions. The improvement parameter, c sw , has 
been chosen so that the leading order discretisation er- 
rors vary quadratically with the lattice spacing |29( . The 
exact Hybrid Monte Carlo (HMC) simulation algorithm 
is used, which avoids finite step-size errors. 

We measure Q using the method of ten cooling 
sweeps are applied using the Wilson gauge action. The 
cooling action has relatively little effect on the topolog- 
ical susceptibility Ten cools strikes a good balance 
between adequate suppression of these ultraviolet dislo- 
cations and excessive destruction of the long range topo- 
logical structure |8j, ll7| . The choice of operator does not 
matter in the continuum limit, at least for the quenched 
theory |30|. and a reflection-symmetrised "twisted pla- 
quette" lattice topological charge operator is used. In 
Fig. n we show the variation of the topological charge as a 
function of Monte Carlo simulation time for three ensem- 
bles over a range of sea quark mass. The histograms of Q 
(using unit width bins, although in the rest of the anal- 
ysis the topological charge is not rounded to the nearest 
integer value) show good agreement with the expected 
Gaussian form. Although (Q) = within statistical er- 
rors, we opt (for consistency with earlier studies) to sub- 
tract terms in (Q) from x- Both this and the choice of 
whether to round Q to integer values affect x by much 
less than one standard deviation. 

This algorithm for measuring x has already been stud- 
ied in the quenched theory (the Wilson gauge action), 
and continuum extrapolation yields C30 = 0.065 (3) and 
c 3 i = -0.28 (4) 0. 

The decorrelation of Q is good (see Fig.^ and |sll^.l25]| 
for details of some other ensembles). More quantitatively, 
we measure the integrated autocorrelation, 7i n t(<3), us- 
ing an implementation of the sliding window method 



|31j , with results shown in Table [H] (scaled to units of 
HMC trajectories) |3!|. With the exceptions of ensembles 
98,12,15 an( A uio, there is weak evidence for a slow increase 
in Ti n t from O(10) HMC trajectories for the least chiral 
ensembles to O(30) trajectories for the most chiral. We 
plot these data for the 16 3 32 ensembles in Fig. |2 (omit- 
ting the least chiral point). The trend is no stronger 
for any particular gauge coupling, /?, but the statisti- 
cal errors and available range of (rom,) 2 make definitive 
statements difficult. Of the outliers, we note that ensem- 
bles 98,15 have only a limited number of trajectories and, 
as we discuss later, q±2 is the smallest physical volume 
studied. 

Barring these possible exceptions, we have confidence 
that the topological susceptibility may be estimated free 
from autocorrelation effects and the statistical errors 
quoted are from a jack-knife analysis using 10 bins. No 
increase in the statistical error estimates for x was found 
with larger bins, as expected when comparing T- mt to the 
bin size. Results for the topological susceptibility are 
given in Tableland plotted against quark mass in Fig. [3] 



A. Discretisation effects 

There is a significant suppression of the topological 
susceptibility in the chiral limit. In Fig. we show the 
leading order chiral expectation of Eqn. (J2J, using the 
physical value for f^. There is approximate agreement 
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FIG. 1: Partial time histories for the topological charge for 
ensembles uio,4 and qio, with normalised histograms of the 
topological charge. 
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for (rom^) 2 < 3, although we note that discretisation 
effects tend to increase x an d exaggerate the range of 
agreement. We do not see clear evidence for \ deviating 
from a rapid decrease with m q at finite lattice spacing as 
in |l7|. although critical slowing prevents exploration of 
sea quarks comparably light with those simulated in that 
study. 

To quantify the effects of lattice artifacts, it is useful 
to consider \ over the full range of sea quark masses. 
Assuming that Eqn. J7J holds, we can use it to perform 
a chiral extrapolation of x whilst allowing for (leading 
order) discretisation effects. We exclude ensemble uio, 
whose statistical error is probably unreliable, and ensem- 
ble q\2 which has a small volume (see on). Ensembles 
58,15 also have large autocorrelation, and we remark that 
excluding them from the fitted data had no statistically 
significant effect. We fit only the dynamical results, with 
the quenched parameters fixed. In principle we could 
avoid this and fit the quenched data as well, assigning all 
such points an arbitrarily large value of (rom^) 2 . Such 
fits do not converge well, however. This was probably 
as a result of the data being arranged as two distinct 
clusters in (rom^) 2 . 

We show the fits to Eqn. in Table Hn with the 
k\ = results (Fit 1) plotted in Fig.[3]for a representative 
sample of lattice spacings, and in Fig. 

It is interesting to note that even for the heaviest sea 
quarks simulated, the discretisation effects act to increase 
X- To see a crossover to the quenched behaviour (where a 
decrease is seen), the fits suggest we need (r m w ) 2 > 25. 
It is thus clear that for all the ensembles the vacuum 
has been qualitatively altered from that of the quenched 
theory by the presence of the sea quarks. As expected, 
topological degrees of freedom appear to be more sensi- 
tive to such changes than most hadronic observables. 



TABLE II: The integrated autocorrelation estimates for the 
topological charge, and topological susceptibility. 



label 


Tint 


4 

r X 




qi 


11 


(2) 


0.0494 


(32) 


U2 


16 


(3) 


0.0445 


(30) 


U3 


30 


(13) 


0.0310 


(38) 


"4 


26 


(8) 


0.0248 


(28) 


U 5 


29 


(15) 


0.0230 


(25) 


96 


12 


(1) 


0.0451 


(32) 


u- 


19 


(4) 


0.0392 


(28) 


qs 


92 


(16) 


0.0302 


(68) 


Ug 


14 


(4) 


0.0433 


(36) 




125 


(40) 


0.0683 


(68) 


fill 


20 


(6) 


0.0369 


(42) 


qi-2 


78 


(35) 


0.0306 


(52) 


qi3 


22 


(8) 


0.0276 


(33) 


qu 


25 


(15) 


0.0356 


(57) 




150 


(22) 


0.0325 


(87) 



TABLE III: Fits to Eqn. Q. Numbers in italic face were 
fixed in the fit. Fit 2 was carried out with two different sets 
of initial values, as described in the text. 



parameter Fit 1 



Fit 2 



coo 
coi 

C30 
C31 

A-i 

2 



xVd-o.f. 



0.0135 (26) 
0.2595 (229) 
0.065 
-0.28 


1.043 



0.0036 (30) 
0.1074 (250) 
0.065 
-0.28 
0.0149 (51) 
0.980 



0.0283 (36) 
1.453 (522) 
0.065 
-0.28 
-0.028 (9) 
0.721 



For ro = 5.0±0.2 (the range studied in |&|9j), these fits 
suggest Co = 0.024 (3)(1)(1), where the errors arise from 
those on coo, coi and the above range of ro respectively. 
This figure agrees very closely with Co = 0.22 (7) obtained 
in the previous work for Eqn. |4"|. and gives confidence 
in the consistency of our fits. 

In addition to the errors arising from the statistical 
variation in the data, there is also a systematic error 
arising from the uncertainty in the quenched behaviour. 
We attempt to account for this by repeating the fits with 
C30 and C31 varying by one standard deviation in each 
direction: coo varies by ±0.003 under this, and it is clear 
that this is a sub-leading effect that can be ignored. 

The fitted c o = 0.0135 (26) corresponds to f„ = 
91.6 (8.8) MeV. This agrees surprisingly well both 
with the experimental (and 2+1-flavour) QCD measure- 
ment (93 MeV) and measurements by other methods on 
N f = 2 SW actions (~ 97.0 (0.3) MeV in these con- 
ventions, rescaled from [12 )■ The leading order lattice 
spacing correction to this increases f„, as expected. The 
size of this correction is consistent with those obtained 
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▲ P=5.29 



I 
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FIG. 2: Integrated autocorrelation estimates for topological 
charge on the 16 3 32 lattice volumes, labeled by gauge cou- 
pling, p. 



5 



for other methods (see Fig. 35 of [32, for instance, al- 
though the different gauge actions and operators will, of 
course, lead to different lattice artifacts). 

It is useful here to consider whether there is any ev- 
idence for \ having a non-zero chiral limit at finite a. 
In Fit 2 of Table IIIII we show the result lifting the re- 
striction k\ = 0. The fits are extremely unstable; the 
first column is the result of starting with all free param- 
eters being zero, whilst the second takes the results of 
Fit 1 as the initial values. This instability probably arises 
from the most chiral data points being clustered around 
ro = 5. Wc conclude that whilst we cannot rule out a 
non-zero chiral limit, the data appears to be described 
more consistently when such a term is not present. For 
the remainder of the paper we therefore concentrate on 
the results from Fit 1. 



B. Finite volume effects 

In 0| the lattice volumes were sufficiently large that 
the topological charge on subvolumes (of size L A ) varied 
independently; indeed a reduction in statistical error on 
X was achieved by exploiting this self-averaging. We did 
not find the same in the this study, and this prompts us 
to consider more carefully whether the topological sus- 
ceptibility exhibits finite volume contamination. 

In general we expect finite volume effects to lead to a 
reduction in x : a small lattice excludes large instantons. 
In the quenched theory this restriction is negligible once 
L/ro > 2.5 |23. With dynamical sea quarks present, 
small volumes can, in addition, see a restoration of chi- 
ral symmetry with the leading order chiral variation now 



0.08 1 , 1 , 1 , 1 , r 




FIG. 3: The variation of the topological susceptibility. To 
guide the eye, the point size is proportional to r^ 2 . The sin- 
gle fitted curve Eqn. |J7) is shown for the continuum limit and 
lattice spacings representative of the plotted data. The expec- 
tation from continuum chiral perturbation theory at leading 
order is also shown. 



being x ^ (mq) Nf °c (m n ) 2Ni . Avoiding this requires 
2Ts 3* 1, although the precise limit is not known. What 
is also not clear is how the behaviour interpolates be- 
tween this and Eqn. J3J as we vary ^ls! we do not address 
this in this study. 

Finally, we need the lattice to be large in units of the 
pion correlation length, Lm^ 3> 1. Spectroscopy suggests 
Lm-x > 5.7 is necessary to render the low- lying hadron 
states free of finite volume effects [34j|. This limit is, 
however, associated more with pion exchange than phase 
structure and is probably unduly stringent for this study. 

In Table|Uwe show all three measures of FVEs. All the 
above limits seem well satisfied and we do not expect, or 
see, any sign of chiral symmetry restoration or significant 
finite volume effects. Certainly we do not see the smaller 
volumes lying consistently below the fitted curve. 

The ensembles 912,13,14 vary only in volume, and no 
strong suppression of x is seen as V is reduced. Nonethe- 
less, we choose to exclude 912, the smallest lattice, from 
the scaling analysis as the volume is particularly small. 

We can attempt to quantify the effects of finite volume 
by repeating the fit with k\ = 0, progressively excluding 
the ensembles with the smallest xls ( u 5, u^...) until the 
fits become unreliable. The fit parameters remained con- 
sistent within statistical errors. We conclude that finite 
volume effects are not significant at this level of statistical 
accuracy if £ls ^10. 

IV. SUMMARY 

We have carried out a systematic study of the topo- 
logical susceptibility, x, m lattice QCD with two light 
quark flavours over a range of lattice spacings and sea 
quark masses using ensembles generated by the UKQCD 
and QCDSF collaborations. We find compelling evidence 
for suppression of x with decreasing quark mass. Finite 




FIG. 4: Combined chiral and continuum extrapolation of the 
topological susceptibility. Solid (blue) points and error bars 
denote data points lying above the fitted surface. Open (red) 
points and dashed error bars denote data points lying below. 
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volume contamination of the results was not seen, and 
x ls ^ 10 is necessary to avoid such effects on the topo- 
logical degrees of freedom. 

We observe a slight (and, in truth, statistically insignif- 
icant) increase in the integrated autocorrelation of the 
topological charge, Q, as the sea quark mass is reduced, 
although this effect is not marked despite a critical slow- 
ing of the HMC algorithm |25j, |2(j . The autocorrelation 
times are, in general, small enough not to bias the jack- 
knife estimates of the statistical errors on our data. 

There have been suggestions that at (5 ~ 5.2 the 
Nf = 2 O(a)-improved SW action exhibits evidence for 
a nearby phase transition 0, l35l |3(j . This transition af- 
fects most strongly observables with the quantum num- 
bers of the vacuum, and we do not expect (or find) a 
strong effect on the topological susceptibility. It may, 
however, induce additional mass dependence in r [36| . 
although we see no signal in r\\ (whose fourth power 
would presumably magnify any such effect, distorting the 
data away from Eqn. (7J). It may also lead to large au- 
tocorrelation (3(|. Whilst a large integrated autocorre- 
lation is seen for one ensemble at large sea quark mass, 
this is at j3 = 5.29 which does not immediately fit the 
hypothesis of |36J. No unexpected trends were seen for 

The chiral variation of the data was compared to a 
theoretically motivated ansatz for the behaviour of \ 
over the full range of sea quark masses. This was ex- 
tended to allow for variation of x with the lattice spacing, 
and a joint chiral/continuum extrapolation was carried 
out for 0.08 < a < 0.12 fm (over a factor of 2 in a 2 ). 
The behaviour was as expected: the continuum limit of 



fir = 91.6 (8.8) MeV was remarkably consistent with 
physical expectations, with discretisation effects tend- 
infcr to increase its value (as seen in other determinations 

The discretisation effects tended to increase the topo- 
logical susceptibility on all the ensembles. This agrees 
with expectations for the chiral limit and is in contrast 
to the quenched theory, where they suppress \- This 
suggests that, relative to the quenched theory, there are 
fundamental differences in the vacuum due to sea quarks, 
even for the largest m q . The topological susceptibility 
is more sensitive to these differences than, for example, 
most quantities in the light hadron spectrum. 

In summary, the topological susceptibility is a sensitive 
probe of the vacuum and chiral properties of the lattice 
action. Although the O(a)-improved SW action formally 
breaks the chiral symmetry at finite lattice spacing, the 
behaviour of \ suggests that the effects arc not significant 
at the lattice spacings simulated. The chiral behaviour of 
this action is comparably good to that of other actions 
currently used in large-scale simulations, and makes it 
a suitable laboratory for studying topologically sensitive 
states, including the rj' meson. 
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